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ABSTRACT 



Aims. Up till now very few dying sources were known, presumably because the dying phase is short at centimeter wavelengths. We 
therefore have tried to improve the statistics on sources that have ceased to be active, or are intermittently active. The latter sources 
would partly consist of a fossil radio plasma left over from an earlier phase of activity, plus a recently restarted core and radio jets. 
Improving the statistics of dying sources will give us a better handle on the evolution of radio sources, in particular the frequency and 
time scales of radio activity. 

Mefhods. We have used the WENSS and NVSS surveys, in order to find sources with steep spectral indices, associated with nearby 
elliptical galaxies. In the cross correlation we presently used only unresolved sources, with flux densities at 1.4 GHz larger than 10 
mjy. The eleven candidates thus obtained were observed with the VLA in various configurations, in order to confirm the steepness of 
the spectra, and to check whether active structures like flat-spectrum cores and jets are present, perhaps at low levels. We estimated 
the duration of the active and relic phases by modelling the integrated radio spectra using the standard models of spectral evolution. 
Results. We have found six dying sources and three restarted sources, while the remaining two candidates remain unresolved also 
with the new VLA data and may be Compact Steep Spectrum sources, with an unusually steep spectrum. The typical age of the active 
phase, as derived by spectral fits, is in the range 10' - 10^ years. For our sample of dying sources, the age of the relic phase is on 
average shorter by an order of magnitude than the active phase. 

Key words. Radio continuum: galaxies - Galaxies: active 



1. Introduction 

Radio sources associated with elliptical galaxies are supplied 
with energy from active galactic nuclei via plasma jets or beams 
(Scheuer 1974; Begelman et al. 1984; Bridle & Perley |1984D . 
In the active stage the total spectra of the radio sources are usu- 
ally well approximated by a power law over a wide range of 
frequencies; spectral breaks at high frequencies, with moderate 
steepening of the spectrum, are also often observed. The break 
frequency is related to the age of the source. 

After this phase, which may last several 10' years, the activ- 
ity in the nuclei stops or falls to such a low level that the plasma 
outflow can no longer be sustained. The radio core, well-defined 
jets and compact hot-spots will disappear, because these struc- 
tures are produced by continuous activity. However, the radio 
lobes may still remain detectable for a long time if they are sub- 
ject only to radiative losses. 

The switch-off of the injection of energetic electrons leads to 
a second high frequency break with an exponential steepening of 
the total spectrum of the radio source (Komissarov & Gubanov 
I1994I I. This second high frequency break is related to the time 
passed since the switch-off of the central engine. Therefore "dy- 
ing" radio galaxies should have steeper high frequency spec- 
tra than normal (active) radio galaxies. Current models of radio 
source physics, based on the assumption that radiative losses are 
dominant, predict the existence of a large population of radio 
sources that are bright at 100 MHz, but have a steep spectrum 
around 1 GHz. However, such sources are not commonly found. 
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According to Giovannini et al. (I 1988I I only few percent of 
objects in a sample of B2 and 3C radio galaxies have the charac- 
teristics of a dying radio galaxy. The first unambiguous example 
of such a source was given by Cordey ( 1987); see also more re- 
cent data on this source in Jamrozy et al. ( 120041 1. Up till now 
only four other sources have been discovered (e.g. Venturi et al. 
1998;Murgiaetal.|2005ll. 

In the literature some possible dying sources, located in rich 
clusters, are mentioned, but these can be interpreted either as 
remnants of dying radio galaxies at the point of being disrupted 
by a particularly turbulent environment (Slee et al. 2001 ), or as 
old radio lobes that are compressed by shock waves (EnBlin & 
Gopal-Krishna [200lT) . However the identification of former host 
galaxies is problematic and thus their association with dying ra- 
dio sources is uncertain. 

It is also possible that radio galaxies may be active intermit- 
tently; in that case one may find fossil radio plasma left over 
from an earlier phase of activity, while newly restarted core 
and radio jets are visible as well. Examples of restarting activ- 
ity in the presence of fossil radio lobes are 3C338 and 3C317 
(Roettiger et al. |1993l Venturi et al. 120041 ); intermittent jet activ- 
ity in 4C29.30 is discussed by Jamrozy et al. ("2007). Restarting 
activity has now also been found in Compact Steep Spectrum 
(CSS) sources (Marecki et al. |2006] l. 

The "double-double radio sources" form an interesting class 
of restarting activity, in which new radio lobes can be seen close 
to the nucleus, while an older pair of lobes is still present far- 
ther out (Schoenmaker et al. 2000). At present about a dozen of 
such sources are known (see Jamrozy et al. 120071 and references 
therein). 
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Dying radio sources and intermittent radio galaxies are more 
easily detectable at low frequencies: a search for this type of 
object should be performed using surveys at frequencies below 
1 GHz. 

In this paper we describe a search for dying sources and dis- 
cuss the properties of the sources we have found so far 

In Sect. |2] we briefly describe the selection of the sample 
and the radio and optical observations. In Sect. [3] we show radio 
contour maps and integrated spectra of the sources in the sample. 
In Sect.|4]we give some comments on individual sources and in 
Sect.|5]we discuss the properties of the dying sources. Finally in 
Sect. |6] we summarize the results obtained in this paper 

All intrinsic parameters (radio power, absolute magnitudes, 
sizes) were calculated with - 71 km s"' Mpc"', Q,„ - 0.3, 

- 0.7. Spectral indices a are defined according to 5 oc v 



2. Selection of the sample and observations 

2.1. The sample 

Since dying radio galaxies and intermittent radio sources 
are more easily detected at low frequencies, the Westerbork 
Northern Sky Survey (WENSS; Rengelink et al. [T997] l at 325 
MHz is particularly well-suited for this purpose. As a starting 
point we have used the work done by De Breuck et al. (20001); 
while their aim was to find high-redshift radio galaxies, we are 
instead interested in steep spectrum radio sources identified with 
bright galaxies: these are our prime dying source candidates. De 
Breu ck et a l. ( |2000l ) correlated the WENSS and NVSS (Condon 
et al. 1 19981 1 catalogues using a search radius of 10 arcsec cen- 
tred on the WENSS positions. They obtained spectral indices 
for 143000 sources with S325MHZ ^ 30 mJy. Only unresolved 
sources (angular cutoff ~ 1 arcmin) were considered. A further 
condition was 5 1400 > 10 mJy. We note that especially this last 
criterion may lead to the loss of many steep spectrum sources, 
because De Breuck et al. (2000) effectively used only that part 
of WENSS with flux densities above ~ 100 mJy. Thus they ob- 
tained a sample of 343 sources with spectral indices Oq^t, > 1-3. 
We plan to extend in the near future the search for dying radio 
sources also to 1400 MHz flux densities below 10 mJy. Optical 
identifications were done using the digitized POSS-I and were 
listed in their table A.4. We checked all the identifications listed 
and estimated galaxy magnitudes following the procedure de- 
scribed in de Ruiter et al. 0998 ). Selecting from list A. 4 only 
sources associated with an elliptical galaxy with nir < 17 we 
obtained 1 1 candidate dying sources. 

2.2. Radio observations 

In order to determine whether the integrated radio spectra of 
these dying source candidates are really steep also at high fre- 
quencies we have observed the sample in snap-shot mode with 
the VLA D-array at 4.8 and 8.4 GHz. The sources are sufficiently 
compact (LAS < 1 arcmin) for being only slightly resolved at 
high frequency and therefore we are able to measure the total 
flux density. However, we need to know whether active radio 
structures such as radio cores, jets and hot spots, are present. In 
order to fill this gap in our knowledge we also performed VLA 
B-array observations, at 1.4, 4.8 and 8.4 GHz. 

Four of the most promising objects were observed with the 
VLA A-array at 327 MHz. In this way we have obtained spec- 
tral information (between 1 .4 GHz and 327 MHz) on the source 
components. 



A summary of the observations, including the array, fre- 
quency, length of observation, beam size and r.m.s. noise is given 
in Table [T] A bandwidth of 50 MHz was used for the B- and D- 
array, while the observations with the A-array were done in line 
mode. 

Calibration and imaging was done using the NRAO 
Astronomical Image Processing System (AIPS) in the standard 
way. 

2.3. Optical observations 

Optical spectroscopy was done with the DOLORES spectro- 
graph installed at the Nasmyth B focus of the 3.5-m Galileo 
telescope on the Roque de los Muchachos in La Palma, Spain. 
The observations were performed in service mode during several 
nights between Jan 29 2004 and Aug 5 2005. We used the LR- 
B grism and a slit size between 1.0 and 1.5 arcsec, depending 
on the seeing conditions during the observations. The dispersion 
was 2.8 A per pixel and the spectral resolution about 1 1 A. The 
typical useful spectral range is ~ 3500 to ~ 8000 A with increas- 
ing fringing beyond 7200 A. The slit orientation was adjusted to 
maximise the number of spectra obtained with a single setup. 

The reduction of the spectra was carried out using the 
LONGSLIT package of NOAO's IRAF reduction softwai-e. A 
bias frame was constructed by averaging 'zero second' expo- 
sures taken at the end of the night. This was subtracted from 
every non-bias frame. The pixel-to-pixel variations were cali- 
brated using internal flats. The sky contribution was removed 
by subtracting a sky spectrum obtained by fitting a polynomial 
to the intensities measured along the spatial direction excluding 
the vicinity of objects. One-dimensional spectra were extracted 
by averaging in the spatial direction over an aperture as large 
as the spatial extent of the source. Wavelength calibration was 
carried out by measuring the positions on the frames of known 
emission lines from a He-lamp and a flux calibration was ob- 
tained from a comparison with spectra of spectro-photometric 
standard stars. 

As the main goal of our spectroscopic observations was the 
determination of the redshifts, the exposure times were kept as 
short as possible. Almost all spectra show the 4000 A break, 
the Ca H and K lines, as well as the Mgl and Nal absorp- 
tion lines. Most sources show no obvious emission lines, except 
WNB 1408. 8 +4630 14, which may have Ha in emission; how- 
ever that part of the spectrum is strongly affected by fringing. 
The resulting redshifts are compiled in Table |2] In addition we 
report the redshifts of some of the ambient objects which were 
covered by the slit (see Table[3]). 

3. Results 

3.1. Radio maps and integrated spectra 

Images of the sources are given in Figs.[T]to[TT] The contours of 
the B-array observations at 1 .4 GHz are superposed on DSS2 red 
images. Contour maps of the A-array observations at 327 MHz, 
shown in Figs. [T2] to [15] are superposed on grey scale images 
of the spectral index (between 327 MHz and 1.4 GHz). D-array 
images are not shown, since all sources are unresolved. 

Measured source, and if possible also component, parame- 
ters are given in Table |2] We give the position of the associated 
optical galaxy and of the source radio core when present. We 
also give the total flux of the source or that of the components 
together with the corresponding observing frequency. Usually 
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Fig. 1. WNB0722. 4+4129. Left: radio contours at 1.4 GHz overlayedon a DSS2 optical image. Contour levels start at 126 /iJy/beam 
(3cr) and scale by V^. The beam is 4'.'3x3'.'7. Right; integrated radio spectrum. The solid line represents the best fit of the synchrotron 
model described in the text. 
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Fig. 2. WNB095 1 .9+6037. Left: radio contours at 1 .4 GHz overlayed on a DSS2 red image. Contour levels start at 78 /vJy/beam (3cr) 
and scale by V2. The beam is 4'.'8 x 3'.'6. Right: integrated radio spectrum. The solid line represents the best fit of the synchrotron 
model described in the text. 



only the central component is detected at high frequency. The 
size of the source corresponds to the full-width-half-maximum 
(FWHM) of a Gaussian fit or to the largest angular size (LAS), 
using the lowest reliable contour 

A dying source should show exponential steepening of its 
total spectrum. The sources of our sample were selected on 
the basis of their spectral steepness between 1.4 GHz and 325 
MHz, but we also collected other spectral information: in addi- 
tion to the D-array observations done at 5 and 8.4 GHz, we also 
made use of CATS (the on-line Astrophysical CATalogs sup- 



port System, at 'http://cats.sao.ru/), thus recovering data at dif- 
ferent frequencies. The 74 MHz flux densities were determined 
directly from NRAO VLSS images that are available on-line (at 
URL: http://lwa.nrl.navy.mil/VLSS/). The results of this search 
are shown in Table |4] where all the available data for each indi- 
vidual source are listed. 

The integrated spectra are shown in the right-hand panels of 
Figs.[I]to[II] 

The total spectra of most of the radio sources in our sam- 
ple show a quasi-exponential cutoff in the frequency range be- 
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Fig. 3. WNB 1023.3+2959. Left: radio contours at 1 .4 GHz overlayed on a DSS2 red image. Contour levels start at 72 juJy/beam (3cr) 
and scale by V2. The beam is 4'.'4 x 3'.'8. Right; integrated radio spectrum. The solid line represents the best fit of the synchrotron 
model described in the text. 
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Fig. 4. WNB 1 127.5+4927. Left: radio contours at 1 .4 GHz overlayed on a DSS2 red image. Contour levels start at 87 /vJy/beam (3cr) 
and scale by V2. The beam is 4'.'5 x 3'.'7. Right: integrated radio spectrum. The solid line represents the best fit of the synchrotron 
model described in the text. 



tween 74 MHz and 8.4 GHz. For some sources the spectral cutoff 
is particularly strong (WNB095 1.9+6037, WNBl 127.5+4927, 
and WNBl 150.0+3749). In some others (WNB 1328. 3 +6520 
and WNB 1458.0+4959) the curvature is less pronounced and the 
integrated spectrum is almost a power law in the considered fre- 
quency range. Finally, WNB1408. 8+4630, WNB1438.0+3720 
and WNB23 17.4+4234 show evident flattening of the integrated 
spectrum at high frequency. This flattening is due to the pres- 
ence of a core-jet component that dominates the spectrum at the 



highest frequencies. This is confirmed by the images at higher 
resolution (not shown), in which only the core is visible. 

3.2. Spectral indices of source components 

WNB095 1.9 +6037, WNBl 127.5+4927, WNBl 150.0+3749, 
and WNB23 17.4+4234 were also observed at 327 MHz with the 
VLA A-array. Spectral index maps, using both 327 MHz and 1 .4 
GHz data at 6 arcsec resolution, are shown in the left-hand pan- 
els of Figs. [T2]to[T5] In the right-hand panels the spectral index 
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Fig. 5. WNBl 150.0+3749. Left: radio contours at 1.4 GHz overlayed on a DSS2 red image. Contour levels start at 102 /iJy/beam 
(3cr) and scale by V^. The beam is 4'.'6x3'.'9. Right; integrated radio spectrum. The solid line represents the best fit of the synchrotron 
model described in the text. 
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Fig. 6. WNB 1257.4+3 137. Left: radio contours at 1.4 GHz overlayed on a DSS2 red image. Contour levels start at 84 /vJy/beam (3cr) 
and scale by V2. The beam is 4'.'6 x 4'.' 1 . Right: integrated radio spectrum. The solid line represents the best fit of the synchrotron 
model described in the text. 



distributions are shown, as computed in boxes along the ridge 
line, or all over the source, in the way shown in the inset. 

In WNB095 1.9+6037 the northern lobe (or tail) has a' = 
1 .2 in the centre, and steepens going outwards in both directions 
to values of q-q^j = 1.9 and Oq'^^, ~ ^-^ '■^^ respective ends. 
The south lobe has a uniform spectral index of Cy^j ~ 1 .7. Note 
that also the central component of this source has a steep spec- 
trum, with a ~ 0.9. In WNBl 127.5+4927 the spectral indices 
of the two lobes are almost constant with values of the order of 
1.3. At most there might be a very slight steepening of no more 



than ~ 0.2 (see Fig. [IS. In WNBl 150.0+3749 the spectral in- 
dex varies from (Xq^-j^ - 1.6 in the source centre to cKq^j = 2.2 in 
the outer parts. In WNB23 17.4+4234 only the head of the tail is 
visible at 327 MHz; it has a},\^ = 0.6. 



4. Comments on individual sources 

- WNB0722.4+4129: a small asymmetric double without a 
core. It is completely resolved out in the B-array observa- 
tions at 5 and 8.4 GHz. 
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Fig. 7. WNB1328. 3+6520. Left: radio contours at 1.4 GHz overlayed on a DSS2 red image. Contour levels start at 180 /iJy/beam 
(3cr) and scale by V^. The beam is 5'.'3 x3'.'6. Right; integrated radio spectrum. The solid line represents the best fit of the synchrotron 
model described in the text. 
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Fig. 8. WNB1408. 8+4630. Left: radio contours at 1.4 GHz overlayed on a DSS2 red image. Contour levels stait at 150 yuJy/beam 
(3a-) and scale by V2. The beam is 4'.'4x3'.'7. Right: integrated radio spectrum. The solid line represents the best fit of the synchrotron 
model described in the text. 



- WNB095 1.9+6037: the central component is detected in the 
B-array observation at 8.4 GHz but not at 4.8 GHz due to the 
higher rms noise on the 4.8 GHz map. Since it has a'('^3 =0.9 
and a^-^ = 1 .2 it is difficult to consider this an active core. At 
327 MHz and 1.4 GHz the source structure is very similar, 
the northern component being slightly more elongated at 327 
MHz. 

- WNB 1023.3+2959: a small asymmetric double without 
core. It is completely resolved out in the B-array observa- 
tions at 5 and 8.4 GHz. 



- WNB 1 127.5+4927: a classical double source with a com- 
pact (apparently background) source superimposed on the 
western lobe. No core has been detected. At 5 and 8.4 GHz 
(B-array) only the background source is detected. The dou- 
ble source structure is very similar at 327 MHz and 1 .4 GHz, 
but the background source is not detected at 327 MHz due to 
its flat spectral index (x; 0.2). 

- WNB 1 150.0+3749: a source with amorphous structure. It 
may be a small double source without a core and is com- 
pletely resolved out in the B-array observations at 5 and 8.4 
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Fig. 9. WNB 1438.0+3720. Left: radio contours at 1 .4 GHz overlayed on a DSS2 red image. Contour levels start at 96 juJy/beam (3cr) 
and scale by V2. The beam is 4'.'4 x 3'.'7. Right; integrated radio spectrum. The solid line represents the best fit of the synchrotron 
model described in the text. 
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Fig. 10. WNB 1458.0+4959. Left: radio contours at 1.4 GHz overlayed on a DSS2 red image. Contour levels start at 150 yuJy/beam 
(3a-) and scale by V2. The beam is 4'.'4x3'.'7. Right: integrated radio spectrum. The solid line represents the best fit of the synchrotron 
model described in the text. 



GHz. At 327 MHz the radio structure is very similar to the 
one at 1.4 GHz. 

- WNB1257.4+3137: a source with amorphous structure at 1.4 
GHz. At 4.8 GHz only the central component is still visible, 
which, however, is completely resolved out at 8.4 GHz. 

- WNB 1328.3+6520: a compact radio source with a very steep 
spectrum. It is detected in all B-array observations. This 
source resembles a Compact Steep Spectrum (CSS) source, 
but with an unusually steep spectrum. 



- WNB 1408.8+4630: apparently this source has one-sided 
morphology. The galaxy position coincides with the bright 
compact component. Much fainter emission extends to the 
south-west. Perhaps this is a core-jet radio source or, alter- 
natively, a small head-tail source. The bright component has 
Q = 0.2 and is probably the radio core. The very steep 
spectrum low frequency emission probably comes from an 
extended component that is not seen at frequencies > 1 GHz. 
However, the integrated spectrum shows flattening at high 
frequencies due to the presence of the flat-spectrum core. 
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Fig. 11. WNB23 17.4+4234. Left: radio contours at 1.4 GHz overlayed on a DSS2 red image. Contour levels start at 96 /iJy/beam 
(3cr) and scale by V^. The beam is 4'.'3x3'.'7. Right; integrated radio spectrum. The solid line represents the best fit of the synchrotron 
model described in the text. 
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Fig. 12. WNB095 1 .9+6037. Left panel: spectral index in grey scale, superposed on radio contours at 327 MHz. Radio contours start 
at 3cr and increase by V2. Right panel: profile of the spectral index along the lobes (see also the text) 



- WNB 1438.0+372: a small bent radio source at 1.4 GHz. A 
core (a^ y = 0.4) is detected at higher frequencies. This radio 
source may be a small wide angle tail (WAT). 

- WNB 1458.0+4959: a very compact radio source at 1.4 and 
5.0 GHz. It is not detected at 8.4 GHz. This source resembles 
a CSS source, but with an unusually steep spectrum. 

- WNB23 17.4+4234: a head-tail source. The tail is very weak 
and only the head is detected at 5 and 8.4 GHz. At 327 MHz 



the tail cannot be seen because the rms noise of the map is 
rather high. The head has Oq^j = 0.7 and ck^q = 0.5. The 
associated galaxy is NGC 7618. 
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Fig. 13. WNB 1 127.5+4927. Left panel: spectral index in grey scale, superposed on radio contours at 327 MHz. Radio contours start 
at 3cr and increase by V2. Right panel: profile of the spectral index along the lobes (see also the text) 
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Fig. 14. WNB 1 150.0+3749. Left panel: spectral index in grey scale, superposed on radio contours at 327 MHz. Radio contours start 
at 3cr and increase by V2. Right panel: profile of the spectral index along the lobes (see also the text) 



5. Discussion 

5.1. Morphological properties and radio cores 

^From a morphological point of view the sources discussed here 
form a rather mixed bag. The only common characteristic is their 



small angular size, but of course this is due to the selection bias 
(< 1 arcmin). Since the sources cover a relatively large redshift 
range (from 0.02 up to about 0.26) the linear sizes show more 
variety, although in a power-linear size plot these sources mix 
with, for example, the general sample of B2 radio galaxies. It 
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Fig. 15. WNB2317.4+4234. Spectral index in grey scale, super- 
posed on radio contours at 327 MHz. Radio contours start at 3cr 
and increase by V2. 



should be noted however, that the morphologies found in this 
sample are very different from those of the B2 radio galaxies, 
but, again, this may be due to the size selection imposed on the 
sample of candidate dying sources. 

A flat spectrum radio core is unambiguously present only in 
a minority of cases: using the higher frequency data we detected 
a radio core in three sources, while we can put upper limits on 
the core luminosity for the other ones. Radio cores might still 
be present at lower levels without flattening the spectrum signif- 
icantly. We find that nine sources have a ratio S foK^'/S < 
7 X 10""^ (of which one is a detection at the level of ~ 5 x 10"^), 
while only two are above this (of order 10"^). ^From the rela- 
tion between core and total luminosity given by Giovannini et 
al. ( 1 19881 1 we would have expected a ratio > 7 x 10"^ for about 
half of the sources. Therefore, even allowing for the selection 
bias introduced by the steepness of the spectra, the cores in the 
present sample are quite weak and this fact strengthens the suspi- 
cion that the majority in our sample may be genuine dying radio 
sources. 



5.2. Equipartition parameters 

We assume that the radio sources contain relativistic particles 
and magnetic fields uniformly distributed and in energy equipar- 
tition conditions. The equipartition parameters (total energy Ueq, 
energy density u^g and magnetic field B^g) are generally com- 
puted assuming that the relativistic particle energies are confined 
between a minimum £„„„ and a maximum /Sma v, corresponding 
to the observable radio frequency range, typically 10 MHz - 100 
GHz (see, e.g., Pacholcz vk 1 1 9701 ) . This choice minimizes the 
source energetics required by the observed radiation in the ra- 
dio band. However, a fixed frequency range corresponds to an 
energy range that depends on the source magnetic field, which 
may change from source to source. Furthermore, as pointed 
out by Brunetti et al. (1997 ) a fixed frequency range computa- 
tion would miss the contribution from lower energy electrons. 



since the £■„„„ corresponding to 10 MHz is larger than 200 
MeV for B^g < 30 jiG; this missing contribution can be very 
large. Because of this we have computed, wherever possible, the 
equipartition parameters assuming a fixed low energy cutoff £■„„„ 
- 10 MeV (for this choice see Brunetti et al. 119971 ). Intrinsic pa- 
rameters, including equipartition parameters, are listed in Table 
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For a power law energy distribution N{E) cc E ^ (6 = 2a+l) 
the equipartition magnetic field is given by the following expres- 
sion 



Beg = lF{a)(l + k) 



L(v*)(yT 



V 



7\-2a 



-E 



\-2a\ 



2a- 1 



I/(3+ff) 



(1) 



where a is the spectral index (assumed to be 0.8), V the source 
volume in kpc^*, v* a reference frequency in GHz, L{v*) the 
monochromatic luminosity at v* in W/Hz, k the ratio of proton 
to electron energy (assumed to be one) and F(a) a quantity in 
the range ^ IQ-^'^ - IQ-^^ for 0.5 < or < 1.0. 

As far as the choice of E^ax is concerned (see Sect. l5.3l for 
details about the parameters used) we note that beyond E{vi,2) 
no electrons are present and the energy contribution in the range 
from E{vh) to E„ax is always much less than that below E{vt). 
Furthermore it turns out that in the equipartition formula E^ax ~ 
E(vh) can be neglected. 

The equipartition parameters computed in this way are re- 
ported in Table |5] A comparison with fixed frequency range cal- 
culations shows that the values of our are larger by up to a 
factor two. 



5.3. Modelling the integrated radio spectra 

The radiative energy losses are assumed to be dominant with 
respect to other processes (e.g. adiabatic losses). The pitch an- 
gles of the radiating electrons are assumed to be continually 
isotropized in a time that is shorter than the radiative time scale. 
According to this assumption the synchrotron energy losses are 
statistically the same for all electrons. After its birth the source 
is supposed to be fuelled at a constant rate (i.e. the continuous 
injection phase) by the nuclear activity, for a duration The 
injected particles are assumed to have a power law energy spec- 
trum N{E) cc E^^'"', which will result in a power law radiation 
spectrum with spectral index ai„j - (dinj - l)/2. In this phase 
the source radio spectrum changes as a function of time in a way 
described by the shift of break frequency vt, to ever lower values 
as the time, f s, increases: 



B 



Vb ^ 



(B2 



(2) 



where B and Bjc are the source magnetic field and the inverse 
Compton equivalent magnetic field, respectively. Below and 
above Vh the spectral indices are respectively ai„j and ai„j+0.5. 

At the time td the power supply from the nucleus is switched- 
off'. After that a new phase of duration fopF begins (i.e. the relic 
phase). A new break frequency Vb2 then appears, beyond which 
the radiation spectrum drops exponentially. This second high fre- 
quency break is related to the first by: 



Vbl = Vb 



- 

VOFF 



(3) 



where fj - + foFF is flie total source age (see e.g. 
Komissarov & Gubanov [T^94l Slee et al. l200TT ). 
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Thus, the above synchrotron model (hereafter CIoff) is de- 
scribed by four parameters: 

i) ai„j, the injection spectral index; 

ii) Vh, the lowest break frequency; 

iii) (off/ ts, the relic to total source age ratio; 

iv) norm, the flux normalization. 

In the CIoFF model the magnetic field strength is assumed to 
be uniform within the source. Slee et al. (2001) and Murgia (in 
prep.) propose a more sophisticated model (MJP), which takes 
into account diffusion of synchrotron electrons in a magnetic 
field with a Gaussian distribution of field strengths. Compared 
to the basic CIqff model the MJP model is characterized by two 
more free parameters: 

v) Brms/Bic, the ratio of the rms to the inverse Compton mag- 
netic field; 

vi) Dgff, the electron diffusion efficiency. 

If the rms magnetic field is stronger than the inverse 
Compton magnetic field (Bims > 5Bic) and the electron diffu- 
sion is low (Deft - 0) the spectrum beyond Vh2 is less steep than 
the exponential drop in the CIqff model. The differences in the 
shape of the spectral cutoff are significant only for extremely 
steep spectrum sources, i.e. if Vb < 100 MHz and foFF/f? > 0.2. 
For this reason we fitted the total spectra of all sources with the 
CIqff model, excluding WNBl 150.0-H3749. In this case only the 
fit obtained by the MJP model was definitely better This is not 
so surprising as WNB1150. 0+3749 is the source with the steep- 
est spectrum of all. The model fit is represented by a solid line 
in the right-hand panels of Figs. 1 to 11, where we also show 
the reduced and the best fit parameters. The most difficult pa- 
rameter to constrain is ai„j, since its determination requires flux 
density measurements at frequencies below 74 MHz, our lowest 
frequency. For this reason, we preferred to keep its value fixed 
to Oinj = 0.7 in the fitting procedure. This value is character- 
istic of the spectra of jets and hot spots of many active radio 
galaxies and can be assumed as a plausible injection spectral in- 
dex. If we used a 



"1J 



- 0.5 instead, we obtain similar results, 
although the average is then slightly worse. Beyond the high- 
frequency spectral cutoff in WNB 1 150.0-1-3749 the spectral drop 
is flatter than exponential. Therefore, this spectrum was fitted 
with an MJP model in which we fixed the diffusion efficiency at 
Deff = and the ratio at Brms/^ic - 2.4. The value of Bims/Sic 
was derived by adopting for B^ms the equipartion magnetic field 
strength Z?eq - 12juG. 

Overall the fits are quite good. The worst cases are 
WNB1328.3+6520 and WNB1458.0+4959, which happen to be 
the two most compact sources of our sample. Apart from these 
two objects, we find that Vb is typically in the range 100 MHz 
- 1 GHz, while Iqff Its > 0.1 confirming that these really are 
radio galaxies in which the fuelling of relativistic electrons has 
been in an off-state by a significant amount of time. The spectra 
of WNB 1328. 3 -1-6520 and WNB1458.0+4959 are steep power 
laws. According to the adopted model we see the aged spectrum 
of the continuous injection phase. Thus, these sources may not 
yet be in the dying phase. 

In general the radio cores are very weak, most having 
^coS"V5m ™' < 3 X The sources WNB1408.8-H4630, 

WNB 1438.0-1-3720 and WNB2317.4H-4234 are different cases. 
Their low frequency spectrum is rather steep, but strongly flat- 
tens above 1.4 GHz, due to the presence of a bright core that 
is clearly detected in the 5 and 8.4 GHz B-array images. These 
cores have a flat spectrum characterized by c < 0.5 between 



these two frequencies and are likely still active. In order to ac- 
count for these components in the spectral fit, we added a power 
law to the CIqff model, using the observed spectral index and 
flux density as the normalization. In these three cases we may be 
observing fading lobes (produced by a previous duty cycle), in 
conjunction with restarting activity in the core. 

5.4. Modelling the spectral index profiles 

In three cases (WNB095 1.9+6037, WNB 1 127.5+4927 and 
WNBl 150.0+3749) the profile of the 0.3 to 1.4 GHz spectral 
index can be traced all along the lobes. The traces are shown 
in the right-hand panels of Figs. 12 to 14. Assuming a constant 
source expansion velocity the spectral index at a given distance 
from the core can be related to the synchrotron age of the elec- 
trons at that location. In particular, if the electrons were injected 
inside the core during the active phase, we expect that the break 
frequency scales as: 



Vh{d) = 



V&(c/max) 



[(t//t/max) ■ (1 - toFF/ts) + toFF/ts? 



(4) 



where the distance d ranges from d - at the core, up to 
d - fimax at the edge of the source. Thus, the break frequency 
along the lobes varies from a minimum of Vb, at d - c/max, up to 
a maximum of Vbi, at c/ = 0. The two limiting break frequencies 
Vb and Vb2 are exactly the same as given in Eqs. ^ and Q, 
respectively. 

WNB095 1.9+6037 has a Z-shape and this complicates the 
interpretation of the spectral index profile: it is not possible to 
identify the regions where the electrons were injected during 
the active phase. On the other hand, in WNBl 127.5+4927 and 
WNBl 150.0+3749 the spectral index steepens going from the 
regions near the core outward. Thus, it is plausible that the elec- 
trons were injected near the core and subsequently flowed out- 
ward. We fitted the spectral index trend in WNBl 127.5+4927 
with the CIqff model and obtained a minimum break frequency 
of Vb - (1800 ± 250) MHz with toFF/t., > 0.79. Such values 
are in very good agreement with those found from the fit of the 
global spectrum. In WNBl 150.0+3749 the fit of the MJP model 
to the observed spectral index trend gives Vb - (260 ± 25) MHz 
and foFF/fj - 0.55 ± 0.05. For this source the CIqff model pro- 
duces a significantly worse fit of the spectral trend than the MJP 
model, in good agreement with the fit obtained for the integrated 
spectrum. We note that a source in which the activity had been 
stopped already for a long time is expected to have Vb2 - Vb and 
thus should have a rather constant spectral index distribution. 
Strong spectral index variations along the lobes can be observed 
only if Vb2 » Vb, i.e. just after the injection has been switched 
off. 



5.5. Radiative ages 

Assuming a constant magnetic field and no expansion after the 
switch-off, the total source age can be calculated from the break 
frequency, Vb'. 



f, = 1590 



B' 



.0.5 



(B2 + B2^)[(l+^)y,]0.5 



(5) 



where the synchrotron age f , is in Myr, the magnetic field in 
juG, the break frequency Vb in GHz, while the inverse Compton 
equivalent field is Bic = 3.25(1 + zf'fJ.G. 
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By adopting the equipartition value Beq for the magnetic field 
strength we can thus derive the synchrotron age fj. Finally, from 
the ratio foFp/f j, which is also given by the fit, we can determine 
the absolute durations of the active and relic phases, and fopp. 
The typical ages of the active phase are in the range 10^ - 10** 
yrs. The relic phase is typically shorter by an order of magni- 
tude, but there are cases in which the duration of the relic phase 
is comparable, if not longer, than that of the active phase (e.g. 
WNB1127.5+4927 and WNB1150.0+3749). It is quite hard to 
get a firm estimate on the relative ages of the active and dying 
phases based purely on source statistics, as we would need un- 
biased and complete samples, constructed in the same way, of 
both dying and active sources. However, we do have some gen- 
eral data on the identification of WENSS sources with nearby 
galaxies (de Ruiter et al., in preparation): preliminary identifi- 
cation work gives us a sample of about 350 radio sources, ex- 
tracted from the WENSS survey, if we follow the same criteria 
(except the steepness of the radio spectrum, of course) as used 
for the dying sources: identification with an elliptical galaxy 
with red magnitude < 17, flux density at 327 MHz > 100 mJy, 
and unresolved at 327 MHz. The ratio of dying to active sources 
would then amount to roughly 2-3 % (8/350), a bit less than the 
~ 10 % implied by the synchrotron calculations given above. 
However, the condition that the source be unresolved (used by 
de Breuck et al. 120001 in their construction of a sample of steep 
spectrum sources) provides un unknown but probably significant 
bias against finding dying sources, since these are expected to be 
on average older and more extended than younger, still active, 
sources. Therefore the discrepancy by a factor of 4 between the 
ratio of dying to active sources estimated from source statistics 
as compared to the value obtained from synchrotron ages, is by 
no means worrying. 

All the ages are given in Table |6] together with the best fit 
parameters. 

It is interesting to compare the radiative ages of the dying 
sources with those of radio galaxies with similar luminosities in 
the low frequency range, that are still in the continuous injection 
phase. Data on active radio galaxies were taken from the sample 
of B2 radio galaxies (Parma et al. 119991 ). For these sources we 
have recomputed the equipartition and the radiative ages follow- 
ing the procedure of Sect. 15.21 The dying sources indeed tend to 
be older than the B2 radio galaxies. In fact, the median age (f,) 
of the dying sources is 63 Myr, against 16 Myr of the B2 radio 
galaxies; all of the dying sources are older than 16 Myr, the me- 
dian age of the B2 radio galaxies. Even if we consider instead the 
median duration of the active phase of the dying sources in our 
sample, we find that this is 48 Myr, three times the age of the B2 
sources. These results are consistent with the scenario in which 
the sources dubbed "dying" indeed represent the last stage in the 
lives of radio sources. 

5.6. The environment of the dying sources 

We investigated the gaseous environment of the dying radio 
galaxies by checking for diffuse X-ray emission in the Rosat All 
Sky Survey (RASS). We extracted the 0.1 - 2.4 keV image of 
a 10' X 10' field around each radio source. The X-ray images 
were corrected for the background and smoothed with a cr = 45" 
Gaussian kernel. 

The RASS count-rate contours are shown in Figs. [T6] and [TtI 
overlayed on the 1.4 GHz VLA B-array images. 

Cross-correlation of the NED database with the RASS X-ray 
images gives the following results: 



- WNB0722.4H-4129. The radio source is located near the cen- 
tre of the Abell cluster A 580, which has also been detected 
in the RASS. 

- WNB095 1.9-1-6037. The RASS image suggests the presence 
of a faint X-ray source coinciding with the dying radio 
source, but no clusters are known within 10' from the radio 
galaxy. 

- WNB1023. 3+2959. No X-ray emission is seen down to the 
sensitivity level of the RASS; no cluster is known within a 
10' radius from the radio galaxy. 

- WNBl 127.5 -h4927. No X-ray emission is seen down to the 
sensitivity level of the RASS; no cluster is known within a 
10' radius from the radio galaxy. 

- WNBl 150.0-H3749. No X-ray emission is seen down to 
the sensitivity level of the RASS; NED reports a cluster 
(NSCJl 15245-H373321) at i:9 from the radio galaxy, but the 
photometric redshift of the cluster is z=0. 1456, while the ra- 
dio galaxy has a spectroscopic redshift of z=0.228. 

- WNB1257.4-H3137. A faint X-ray source is present in the 
RASS image at the position of the dying galaxy. NED lists a 
cluster (NSCJ125947+312215) at 1.'7 from the radio source, 
with a photometric redshift of z= 0.0585; this is in agree- 
ment with the spectroscopic redshift of the radio source, 
z=0.05172. 

- WNB1328. 3-1-6520. An extended X-ray source is present in 
the RASS image, a few arcminutes north of the radio source 
position. The nearest cluster, Zw 1328. 9-1-6518, lies4' south- 
east of the radio source. It is therefore not clear whether the 
X-ray emission is associated with this cluster 

- WNB1408.8-H4630. Faint X-ray emission is present in the 
RASS image, at the position of the radio source. The 
X-ray emission is probably associated with the cluster 
NSCJ141045-H461634 (photometric redshift of z=0.178). 
However the redshift of the radio source is 0.13255. 

- WNB1438.0-H3720. The radio source belongs to the cluster 
MACSJ1440-H3707, which has been detected in the RASS. 

- WNB 1458.0-1-4959. The radio source is located in the clus- 
ter Abell A 2011. The RASS image shows a peak of X-ray 
emission in correspondence with the radio source and diffu- 
sion emission about 3' to the east. It is unclear whether the 
source is located at the cluster centre or in the periphery; 

- WNB2317.4-H4234. The RASS image shows strong diffuse 
X-ray emission, confirmed by a Chandra observation (see 
bottom-right panel of Fig. [T3l l. at the position of the radio 
source and coincident with the galaxy NGC 7618. 

No firm conclusions can be drawn about the environment 
of dying sources, due to the small number of objects involved. 
Although about half of the dying sources are located in clus- 
ters, there are some (WNB 1023.3+2959, WNB 1 150.0+3749 
and WNB 1 127.5+4927) which appear to be isolated. 

6. Conclusions and Summary 

The definition of dying source is based on spectral shape and the 
absence (or at least weakness) of a radio core powering large- 
scale jets. Using this definition our search for dying sources 
started with eleven candidates, which in the end revealed six 
truly dying sources, and three more in which a radio core de- 
tected at higher frequencies suggests activity that has recently 
restarted. 

We emphasize that this is the total number of dying sources 
present in the entire WENSS catalogue, at least of those with an- 
gular sizes less than about one arcmin, flux density at 1 .4 GHz 
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> 10 mJy in the NVSS catalogue, and identified with galaxies 
with m,- < 17. Because of these restrictions we may very well 
have missed sources with very steep spectra (i.e. those with a 
20 cm flux below the adopted 10 mJy limit) and also the more 
extended sources. We are constructing a new sample of steep 
spectrum sources, always using WENSS, NVSS and an optical 
magnitude limit of 17 in the red, but this time without any re- 
strictions on angular size restriction and 20 cm flux density (de 
Ruiter et al. in prep.). Until then it is very hard to give an es- 
timate of the frequency of dying sources in general radio cata- 
logues, and in particular a statistical estimate of the relative time 
scales of the dying and active phases, to be compared with the 
foFF / ts directly obtained by the model fits of the radio spectra 
(see Table|6l). 

The typical ages of the active phase derived by the spectral 
fits are in the range 10^ - 10^ yrs. The relic phase is typically 
shorter by an order of magnitude, but there are cases in which 
its duration is comparable, if not longer, than that of the active 
phase (e.g. WNB1127.5H-4927 and WNB 1150.0-1-3749). A dy- 
ing source tends to be older than a general B2 radio galaxy. In 
fact, the median age of the dying sources is 63 Myr, against 16 
Myr of the B2 radio galaxies; all of the dying sources discussed 
here are older than these 16 Myr. As we argued in the previ- 
ous section, it is difficult to get an estimate of the relative ages 
from source statistics, due to the present lack of well defined 
samples of both dying and still active sources. Preliminary data 
suggest that synchrotron ages predict about four times as many 
dying sources than actually found, but if we take into account the 
possible biases in our search procedure, plus the intrinsic uncer- 
tainties in the synchrotron calculations, we do not consider this 
a serious discrepancy. 

Two of the sources are only marginally resolved at our best 
resolution. Their linear sizes are < 10 kpc, so they formally be- 
long to the class of Compact Steep-Spectrum Sources, although 
they are less powerful than the classical members found in the 
3CR and PW catalogues. 

Although about half of the dying sources are located in clus- 
ters, there are some which appear to be isolated. Note that the 
three dying sources found in the WENSS minisurvey sample 
(de Ruiter et al. 1998) are all in a cluster (Murgia et al. 2005]). 
Considering the small numbers involved, no firm conclusions 
can be drawn about the environment of dying sources. 

The radio luminosities at 325 MHz are close to the break in 
the luminosity function. According to the spectral evolution sce- 
nario we have assumed that the radio power before the switch- 
off of the central engine was somewhat higher and has decreased 
due to the radiation losses, but by no more than a factor 2-3. 
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Fig. 16. The VLA 1.4 GHz B-array image (grey scale), superposed on the RASS X-ray contours. Contour levels start from 10 ^ 
cts/s (3cr) and scale by V2. 
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Fig. 17. First five panels: the VLA 1.4 GHz B-array image (grey scale), superposed on the RASS X-ray contours. Contour levels 
start from 10"^ cts/s (3cr) and scale by V2. Panel at bottom-right: overlay of the Chandra image and radio contours at 1.4 GHz 
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Table 2. Observational data 
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DEC 


redshift 
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Table 3. Additional redshifts 



Parma et al.: Dying sources 



Name 


RA(J2000) 

hm s 


DEC(J2000) 

Off/ 


Redshift 


WNB095 1.9+6037 
WNB 1150.5+3749 
WNB 1438.0+3720 
WNB 1458.0+4959 


09 55 31.32 
11 52 35.97 
14 40 03.83 
14 59 43.44 


60 23 15.3 
37 32 51.7 
37 06 55.0 
49 47 35.2 


0.20603 
0.22867 
0.09483 
0.16545 



Table 4. Integrated flux densities 
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Name 


5 74 




5 130 
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5 1.4 


^48 


^84 
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0.12 


2,6 


w iNo L\jz.j.j-Tz.yjy 








o 1 






1 1 n 


U.oU 


u.zo 


A 


WNB 1127.5+4927 


495 


364 


228 


132 




113 


17.6 


0.28 


<0.09 


5,6,8 


WNB 1150.0+3749 


4600 


1582 


620 


414 


341 


250 


17.1 


0.43 


<0.08 


5,6,7,8 


WNB1257.4+3137 


792 


540 




268 






34.0 


3.29 


1.09 


4 


WNB 1328.3+6520 


369 






74 






10.5 


1.72 


1.09 




WNB 1408.8+4630 


527 


467 




103 






14.2 


6.51 


5.30 


5 


WNB 1438.0+3720 
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40.2 


4.25 
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Data at 4.8 GHz and 8.4 GHz are from the present paper; data at 325 MHz and 1.4 GHz are from De Breuck et al. ( |2000t . References for other 
frequencies: 

1- Hales et al. ( 19881: 6CII, 151 MHz 

2- Hales et al. ( 19901: 6C1I1, 151 MHz 

3- Hales et al. ( 1993): 6CV1, 151 MHz 

4- Waldram et al. ( 19961: 7C4, 151 MHz 

5- Riley et al. ( 19991: 7C, 151 MHz 

6- Zhang et al 1 1997 1: Miyun, 232 MHz 

7- Douglas et al . (1996j : TXS, 365 MHz 

8- Ficarra et al. l fT985T l: B3.1, 408 MHz 
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Table 5. Intrinsic parameters 



Parma et al.: Dying sources 



Name 


redshift 


M, 


log(P,/0.33MHz) 


Linear Size 
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0.4 
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25.01 


0.60 
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33.6 
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-23.5 
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0.09661 


-22.9 


24.70 
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1.5 


13.6 


WNB 1458.0+4959 


0.16782 


-23.7 
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6.0 X 3.1 


1.6 


44.8 


WNB23 17.4+4234 


0.01730 


-21.0 


22.99 


11 x7 







Table 6. Spectral parameters 



Parma et al.: Dying sources 



Name 




Vh 

MHz 


tow Its 


MHz 


Myr 


tci 

Myr 


fOFF 

Myr 


Notes 


WNB0722.4+4129 


1.2 
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12+""* 
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WNB 1328.3+6520 
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28 


0.6 
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15 
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Notes: 

(1) WNB 1127.5+4927; Integrated CW model. 

(2) WNB 1127.5+4927; Sp. index profiles, JP model. 

(3) WNB 1150.0+3749; Integrated MJP model with fi/Sic=2.4 fixed (i.e. B = \lnG). 

(4) WNB 1150.0+3749; Sp. index profiles, JP model with B/Bic=2A fixed (i.e. B = 12/iG). 



